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Abstract 

The distortions in the thermal energy spectra for neutrinos produced in a su- 
pernova when a resonant oscillation, MSW effect, occurs are determined. In order 
to show this effect for some relevant and representative examples of unified gauge 
models, we have chosen S'O(IO), and SU{5)susY , SO{10)susY with a particular 
scheme for fermion masses (DHR model). The analysis has been performed for two 
choices of neutrinos parameters, predicted by the above models, and capable to ex- 
plain the solar neutrino problem. In both cases one observes a strong distortion in 
the electron neutrino energy spectrum. This effect, computed for a wide range of 
SO (10) SUSY models has produced the same results of the previous supersymmetric 
ones. 
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1 Introduction 



The long-base neutrinos physics certainly represents a suitable arena to reveal the presence 
of fundamental effects, which due to their weakness would be at the moment out of the 
range of accelerators experiments. In particular, the explanation of the so-called solar 
neutrino problem seems to suggest the presence of non-vanishing neutrino masses, and 
of a corresponding mixing a la Cabibbo-Kobayashi-Maskawa also in the leptonic sector 
of the standard model. Interestingly, since similar predictions are naturally obtained in 
the framework of unified gauge models for fundamental interactions, one can think to use 
the neutrino measurements to test the validity of such models 0, ^. 

In general, if one assumes massive neutrinos, the flavour eigenstates do not coincide 
with the mass eigenstates, and thus the phenomenon of neutrino oscillations can arise 
In fact, neutrinos are produced by weak interactions in flavour eigenstates, whereas 
they propagates in space as energy eigenstates (mass eigenstates), and these two basis do 
not coincide. Hence, there is a non-vanishing probability that starting with neutrinos of 
given flavour, at a certain distance from the source one can detect a neutrino of different 
flavour. Unfortunately, for solar neutrinos, this simple oscillation mechanism (vacuum 
oscillations) is almost independent of the neutrinos energy, and thus it hardly reproduces 
the experimental defect in the emission, which seems to show a quite strong energy 
dependence [|I|. A possible solution for the solar neutrino problem is then provided by 
the so-called resonant oscillation mechanism, the Mikheyev-Smirnov-Wolfenstein effect 
(MSW) §. 

According to the MSW mechanism, when neutrinos propagate in a medium rich of 
electrons, neutrino oscillations can receive a resonant enhancement due to the different 
interactions experienced by i^e's (through both charged and neutral current) with respect 
to z/^'s and z/^-'s (through only neutral current) with the medium itself. 

Let us consider for simplicity only two neutrino species, and with x = /i or r. In 
this case if we denote with E the neutrino energy, the enhanced oscillation occurs only if 
the following resonance condition is satisfied |^ 

^cos(2^) = V2GpN, , (1) 

where Am? is the squared masses difference between the two neutrinos, 9 the vacuum 
mixing angle, N^, the electron number density of the medium, and Gp the Fermi coupling 
constant. 

As one can see from (|l]), the resonance condition if applied to astrophysical frameworks, 
would allow to scan very low values of Am^ and 6 by using quite low energy neutrinos, 
since very high densities can be achieved for large regions. One of these situations is 
certainly provided during a supernova explosion. 

In this paper, we study how the MSW effect modifies the neutrinos energy spectra 
coming out from a supernova. The predictions are obtained for a wide range of unified 
gauge models. In particular, we have analyzed the S'0( 10) predictions, as the most rele- 
vant example of non supersymmetric unified model still compatible with the experimental 
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bounds on proton lifetime 0, ^. Furthermore, for the sake of completeness, we have also 
considered SU{5) and 5*0(10) supersymmetric unified models ^ 0, |]. 

The paper is organized as follows. Section 2 is devoted to the study of the MSW 
effect for medium with varying density. In particular, the relevant quantities are defined 
and the expressions for neutrino spectra are obtained. In Section 3 we review the main 
predictions for the neutrino physics corresponding to the still viable unified gauge models, 
supersymmetric and not, and for these theories we compute the neutrinos spectra. Then 
in section 4 we give our conclusions and remarks. 



2 MSW effect in supernovae 



In media with varying density the probability of neutrino conversion depends on the rate 
of such a variation. Within the adiabatic approximation 0, namely under the assump- 
tion that the medium density crossed by neutrinos changes very slowly, one predicts an 
high rate of conversion. Quantitatively, this approximation is valid when the adiabaticity 
parameter 



7 = 



Am^ sin2(2e) 



2E cos(2^) 



1 



dr 



(2) 



computed at the resonance point, is greater than 1. In the opposite case, the survival 
probability would be of the order of unity, so that no appreciable conversion can be 
obtained. For neutrinos propagating in a given medium it is useful to introduce the 
adiabatic and non-adiabatic threshold energies 
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(4) 



where N^^'""^ is the electron density at the neutrinos production pointQ. The former indi- 
cates the transition between no matter effects and resonant conversion, described by the 
adiabatic approximation, while the latter represents the energy scale for the transition 
between the adiabatic and non-adiabatic regions. 

The density range in a supernova is very large; it goes from ~ 10~^ g/cm^ in the 
external envelope up to ~ 10^^ g/cm^ in the dense core. For our purposes, as a reasonable 
radial density profile one can use the following scaling law |Q |10 



P 



Po 



Rq 



(5) 



with Po — 3.5 X 10^'^g/cm^, and Rq ~ 1.02 x 10'' cm. The quantity Ro defines the so-called 
neutrinosphere, which represents the bounding surface of the region in which neutrinos of 



^For simplicity, one can assume that this point essentiaUy coincides with the neutrinosphere. 
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a given flavour are in thermal equilibrium (eq. is approximatively valid for r > Rq). 
Outside the neutrinosphere, the electron fraction number Yg can be assumed almost con- 
stant with respect to the radius, and thus the value Yg = 0.42 can be used. Since inside 
the neutrinosphere the resonance conditions are not satisfied, for a wide range of uni- 
fied models, we can assume in first approximation that on this surface the neutrinos are 
emitted with an almost Fermi-Dirac distribution, and with the different flavors equally 
populated. However, flavour changing effects occurring inside the neutrinosphere have 



been proposed in literature |Tl| , but in any case they would represent further corrections 
with respect to the bulk phenomenon we want to describe. 

Consequently, from (||) it follows that neutrinos diffusing out the neutrinosphere ex- 
perience two MSW resonances: the first resonance, between z/g and z/,- , at a density of 
about 10^ glcrn? and the second one, between z/g and z/^ , at about 10^ g/crn?. These 
values refer to average neutrino energies i?,^^ ^ 10 MeV , E^^^^^ ~ 20 MeV , and to the 
indicative values Am^^ ~ 10~^eV^^, Am^^ ~ 25 eV^^ for the squared masses differences. 
The evolution of the different neutrino flavors can be simply summarized by the following 
diagrams 

> (6) 



> 
















> T^e 


> 













(7) 



Vr { > i^u > \P) 

> 

Electron-neutrinos may oscillate into z/^ at the first resonance or into z/^ at the second 
resonance, whereas muon-neutrinos can oscillate only into u^. at the second resonance. 
Finally, tau-neutrinos may oscillate into at the first resonance and the latter can further 
oscillate into z/^ at the second resonance. 

To study the evolution of neutrino flavors it is crucial to know how these resonance 
regions are crossed, namely if adiabatically or non-adiabatically. 



2.1 The z/g survival probability and spectra distortion 

In order to obtain the energy spectra of the emitted neutrinos from a supernova, let us 
focus our attention on neutrinos produced in the thermal phase, which diffuse out the 
neutrinosphere with a spectrum approximatively given by |jlO|, p!2| p| 



0.5546 



, E 

1 + exp ( — 



(9) 



Note that we assume neutrino masses negligible with respect to T^. 
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where the index i denotes the particular neutrino specie. Since the production and scat- 
tering cross-sections for electron-neutrinos are larger than for the other flavors, z/g are 
produced a bit more copiously with respect to the other ones. Thus, their neutrinosphere 
is larger than that for , v-r . This implies that the temperature Tj in (|^) for the Ve - 
sphere is lower than the one of z/^ , i/^ . Furthermore, since and Vr are produced and 
scatter on the surrounding matter only through neutral currents, they have identical spec- 
tra. Obviously, since i/^ , , and z/g, J^^i, T^r are produced in pairs, the magnitude of 
neutrino and antineutrino fluxes are equal for each flavour. For the temperature of Vf. and 



, Vt- neutrinosphere we adopt the typical values [|10 
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Te ~ 3MeV , (10) 
= Tr - 6MeV . (11) 

By indicating with and = = n^^ the spectra of z/g and z/^ , Ur at the relative 
neutrinospheres, according to the oscillation scheme previously outlined, we have that the 
spectra of the different neutrino flavors emerging from a supernova are 

>u,) + P{v, ^ z/g )] nl^ , (12) 
. z/^ ) + P{Pr -^v^ )] nl^ , (13) 

^rX , (14) 

where Piva z/^) is the probability for the transition — > Z//3. Note that, from the 
experimental point of view, due to the relative low energy of the neutrinos considered 
(O-i-50 MeV\ one cannot distinguish between z/^ and , thus the only relevant quantities 
are and + riy^ . Consequently, by using the unitarity conditions on the transition 
probabilities (see diagrams (|^)-(|D) we obtain 

n,^ = P(z/g ^ z/g ) nl^ + [1 - P(z/g ^ z/g )] , (15) 

n., + n,^ = [1 - P{ye ^ Ve )] nl^ + [1 + P(y, ^ z/g )] n^^ . (16) 

From the above equations it is evident that the only relevant quantity to be computed from 
the MSW theory is the probability that z/g survives to both z/g — > z/^ and z/g resonances. 
Since for a wide range of unified models the two resonances are well separated we can 
write 

P(z/g -> z/g ) ~ PeAVe z/g ) Pg^(z/e Z/g ) , (17) 

where Periye —>■ i^e ) and Pefi{i^e ) are the survival probabilities for z/g 's to the first 
and second resonance. These quantities can be calculated by using a simple two flavour 
approach. According to this we have |P 

1 /I 



Pe^Ve^Ve) = 2 + U ~ ^eM)cOs(20g^)cOs(2e^ 
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and an analogous form for Per(^'e ^ J^e )• Here 9m denotes the effective mixing angle in 
matter calculated at the neutrinosphere. It is given by 



cos(26'„ 



sign(EA - E) 







l + tan2(2^) 



-1/2 



E^E, 



E = Ea 



(19) 



where P'^ is the probability of level crossing at the resonance. For the density scaling law 
(I), it has the expression 



P'^ = exp 



^ypie) 



(20) 



where the function F{6) is approximatively given by 



F{9) ~ (l - tan\e)) 1 1 + - 



log (1 - tan^iO)) + 1 



1 + tan2(^) 
tan\e) 



log(l + tan2(0) 



and E^f^ is the effective non-adiabatic threshold energy for density scaling as 1/r^ 



(21) 
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3 Supernova neutrino spectra for unified gauge mod- 
els 

Let us apply the previous results to some relevant example of unified gauge models. In 
particular, we will consider the unified gauge model based on SO (10) 0] which among 
the non supersymmetric GUT theories, represents the most appealing unification scheme. 
Furthermore, we will extend our analysis to the SUSY version of SU{5) and SO{10) 
recently extensively discussed in literature @, 0, |H . 



3.1 5*0(10) GUT models and neutrino predictions 

As already noted in the previous section, the MSW effect strongly depends on the neutrino 
mixing angle 6 and on Am^ . For Ue — >■ transitions, however, one can restrict the values 
of these quantities by requiring that MSW explains the solar neutrino problem [|l|. This 



analysis, performed by one of the present authors pSl has given the following results 



Aml^ ~ 6 X 10"^ eV^ , (23) 
sin2(2^e/.) ^ 7 X 10"^ . (24) 
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A value for the mass of z/^ of the order of the square root of the r.h.s. of Eq. (^) (and 
a value for the mass of Vf. several order of magnitude smaller) is found in 5*0(10) models 
with intermediate symmetry 5'?7(4)p5(or 5'f/(3)c ® U{1)b-l) ® SU{2)l ® SU{2)r. In 
general, a mass for z/^ about rn^/im?^ larger than the one of is expected 0]. This 
imphes, provided mp^ « miy^ , 

mu^ ^35eV , (25) 
larger than ~ 5 eV which is the value found if one identifies z/^ as the hot component of 



dark matter |]T5l. However, the value given by Eq. (PSQ is obtained by assuming equal 
Majorana masses for u^r and UrR- Thus, the factor 7 of discrepancy may be just ascribed 
to this ansatz. For this reason one can safely take for rriiy^ the value 5 eV. 

By applying Eq. (H) this value corresponds to an adiabatic threshold for the z/g -v^ 
transition in the matter of about 10~^ MeV (since we expect small mixing angle we take 
cos(2^eT) — !)• The efficiency in the conversion Ur depends on the value of the 

mixing angle between Ue and . To be more precise, it can be large if the threshold of 
the non adiabatic region is larger than the adiabatic threshold , or small in the opposite 
case, if the adiabatic region vanishes. 

Let us denote with O^t the value of the mixing angle for which the two thresholds (^ and 
(p^) coincide. With the density scaling given by (^ we get 

sin2(2^er) = 3.25 x lO'^ . (26) 

Therefore, for > ^er, the MSW conversion z/g Vj. will be important, otherwise it 
will be small. Then in the following analysis we will consider the two values sin^(2^e^) = 
3.25 10~^ and 3.25 x 10~^°, representatives of the two widely different regimes discussed 
above. 

By using the above values for neutrino parameters we can apply the MSW theory to 
obtain the spectra of emitted neutrinos from a supernova. In particular. Figure 1 shows 
the adiabaticity parameter 7 for v^. z/^ (solid line) and Vf. —>■ Vj. transitions (dashed line), 
respectively, versus the neutrino energy. Note that, the dashed line corresponds to the 
larger value for sin^(26'eT-), since to the other choice just would correspond a line two order 
of magnitude lower. 

The plots show clearly that the Ve resonance is crossed much less adiabatically 
than Vf. — i> z/^ if values for sin^(26'er) << 3.25 x 10~^ are considered. Hence in these cases, 
the transition z/g — > v^i is expected to be more effective than v^. ^ . 

In Figures 2 and 3, 4 the v^. survival probabilities, corresponding to the resonant 
transitions z/g and z/g respectively, are plotted versus the neutrino energy. 

From Figure 3 the extreme non-adiabatic nature of the v^. — > v^. transition for sin^(26'e^) ~ 
3.25 X 10"^° is clear. 

In terms of the survival probability one can calculate the modifications of neutrino 
energy spectra as induced by MSW oscillations. In Figures 5 and 6 are plotted ra^J^ , 
Ui,^ and + , n^^ + riy^ respectively as given by equations (plSj), (|1^). In particular, 
the dashed and dotted lines represents the 5*0(10) predictions for the values sin^(26'eT) = 
3.25 X 10"^° and 3.25 x 10"^ respectively. 
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As one can see from these Figures, the Vf. spectrum becomes harder toward higher en- 
ergies (more easily detectable) while , emerge from a supernova with lower energies. 
The integral flux of is also partially reduced of about 16%. 

However, the possible values of the neutrino mixing angles in 5*0(10) models with 
fermion couplings to a 126 and a complex 10 , as it is the case of the models discussed in 
Ref. [0], has been widely discussed by Babu and Mohapatra ||^. They conclude that it 
is possible to get the z/g -v^ z/^ MSW mechanism advocated to explain the solar neutrino 
experiments with sin^(26'er) ~ 10~^, much larger than the r.h.s. of Eq. (p6|). In this case 
the Vf. Vt transition in supernovae is extremely efficient and thus the corresponding 
Periye ) Eq- (Q) almost Vanishing. Thus, the expressions of the neutrino 
spectra, as can be seen from ([l^) -(p!^, take the simple form 

n.^ = nl , (27) 
nu, = nl , (28) 
n,^ = nl . (29) 

It is interesting to observe that in this case the u^- and distributions just interchange, 
whereas the one remains unmodified. This result can be easily understood by observing 
that for such an efficient -v^ i^t transition, all the final Ue come either from the former 
Ur , not subsequently transformed in z/^ , or from the initial z/^ just transformed in z/g . 
Thus, since the initial distributions of z/^ and Ur are assumed equal, this yields to the 
result (|23)-(^. 



3.2 SU{5)susY and SO{10)susy unified models and neutrino pa- 
rameters 

As an example of a supersymmetric grand unified model, we consider the scheme discussed 
by Dimopoulos, Hall and Raby (DHR) [^] which is based on the Georgi-Jarlskog texture 
p!6| for fermion mass matrices at the GUT scale. A peculiarity of the DHR model is that 



it can be realized both in SU{5)susy and in SO{10)susy GUT models (obviously, the 
Higgs multiplets in the two cases are different). The predictions for the neutrino mass 
ratios of this model are 



— ^ = 9 ^4 ^ ^ 3Q 
mv^ rrit 

^ = 4(^3K'.)-^ , (31) 



where k = 1 or —1/3 (depending on particular choices in the DHR framework, see 0), 
and = rjc (V"J,)^ mt/mc, fjc being the QCD renormalization enhancement factor of 
the charm mass between rrit and rric- The three leptonic mixing angles {6[, 62, 6'^) are 
expressed in terms of the three quark mixing angles (^i, 62, 63) by 

sin(^;) = ~ sm{9,) , (32) 
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sin(e^) = — sm(^2) , (33) 
sin(^3) = 2kri3 sm(^3) . (34) 

In particular, 



— + --^J ^ cos(0) , (35) 

Oer ^ sm{0[) sm{e'^) (36) 
where (j) is the pliase of the leptonic CKM matrix. 

Inserting in the above formulas the known experimental parameters, one obtains 



_ J 3100 ± 1000 I 



[ 38± 12 II 

mu^ _ f 208±42 I 
I 1870 ±370 II 



(37) 
(38) 



,.,2.2. ^ _ J (1-7 ±0.2)10-2 I 

sm [ZUe^,) - I (9.0 ±4.3)10-2 II ' ^^^^ 
• 2/0/. X f (1.3 ±0.3)10-6 I 

sm^(2^e.) = I [1.4 ±0.3)10- II ' ^40) 
where I and II refer to the different values of k. 

In this framework, as done in the previous section, one can fix the ^/i mixing angle 
and squared masses difference to those values capable to explain the solar neutrino prob- 
lem via the MSW mechanism (sin2(20eM) = 6x10-^^2x10-2, Am^^ = (4^-9)10"^ eV"^ 
, see [|l]). This requirement selects the model I with the following squared masses differ- 
ences and mixing angles 

Aml^ ~ 6 X 10-6 e1/2 ^ (4I) 

Aml^ - 0.27 eV^ , (42) 

sin2(2^,^) ~ 1.7x10-2 , (43) 

sin2(2^er) ^ 1-3 X 10-6 . (44) 

In Figures 5 and 6 the predictions for the neutrino spectra for the DHR model, which is 
compatible with both SU{5)susy and SO{10)susy, are shown (dashed-dotted lines). 

As a possible alternative scenario, Lee and Mohapatra ^ have also considered the 
possibility of SUSY 50(10), directly broken at a scale ~ lO^^ GeV to SU{3) (g) SU{2) O 
f/(l). In this scheme one expects Majorana masses for the z/^j's at the highest scale and 
it is the u-r the neutrino expected to have a mass around 10"^ eV. In that case the mass 
of Ufj^ is too small to cause a <S=4> z/^ MSW conversion in an experimentally interesting 
range of neutrino energies. In this case, in supernova one should have only the conversion 
z/g -vv- z^T-, responsible for the solution of the solar neutrino problem. The situation would be 
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very similar, as long as the distortions of neutrino spectrum are concerned, to a situation, 
where the MSW z/g ■<=^ z/^ conversion is responsible for the solution of the solar neutrino 
problem and the mixing angle Oer is smaller than Oer- 

4 Conclusions 

In this paper we have analyzed the distortions induced on the thermal energy spectra 
of Ve and + emitted from a supernova due to MSW effect. In order to study this 
phenomenon for some relevant examples of unification scenarios, we consider as possible 
choices of the neutrinos parameters, namely masses and mixings, the ones corresponding 
to the 5*0(10) gauge group 0, |^, and the SU{b)susY^ S0{1Q)susy unified models in 
the scheme DHR 0, respectively. Other different S0{1Q)susy models have been also 
considered |0] , and interestingly, as far as the neutrino spectra predictions are concerned, 
it yields to results almost equal to the previous supersymmetric ones. Note that, in all 
cases Am^^ and sin^(26'e;i) have been chosen to be compatible with the explanation of 
the solar neutrino problem [|l|. 

According to the typical values for the density distribution in a supernova, out of the 
neutrinosphere, the condition of resonance is satisfied for the two flavour transitions z/g — > 
z/^ and z/g — > z/,-, respectively. By computing the corresponding adiabaticity parameter 7, 
one observes that, at least for 5*0(10) GUT theories, the resonance z/g — * Vt is crossed 
less adiabatically than z/g z/^ if one takes values sin^(2^e^) << 3.25 x 10~^, and thus 
the latter results the most efficient. The survival probabilities for the electron- neutrinos 
corresponding to the above resonances are then obtained. 

Finally, the distortion of the energy spectra are obtained for the Ve and z/^ + z/,-, emitted 
from the supernova neutrinospheres in pure thermal equilibrium distributions. Remark- 
ably, for both models 5*0(10) and SU {^)susy, S0{1Q)susy in DHR scheme, the distortion 
in the z/g distribution is quite relevant, and moves neutrinos from the low energy region of 
the spectrum to the high energy sector. As a result of this analysis, one observes that the 
reduction of z/g flux for 5*0(10) and for the DHR model is of the order 16%. Thus it is 
essentially model independent. This feature can be easily understood by observing that 
the relevant (since almost adiabatic) neutrino transition is z/g — > z/^ , whose parameters 
are fixed from the MSW explanation for the solar neutrino problem. 
The distortion in the z/g distribution becomes extreme when the choice of neutrino pa- 
rameteres discussed in the papers of Babu and Mohapatra |^ is taken. In this case, in 
fact, since the survival probability Per(z^g — * z/g) is almost vanishing, we have the exact 
interchange of z/g and v^. distributions. 
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Fig. 1 




E (MeV) 



Figure 1: The adiabaticity parameter 7 for t^ii (solid line) and for i/g — transitions 
(dashed line) versus the neutrino energy is reported, for the SO{10) unified model |0. 
Note that for the 5*0(10) case the line corresponds to sin^(26'er) = 3.25 10~^. 
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Fig. 2 



^ 0.8 
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Figure 2: The survival probability Pe^{i^e ~^ ^e) as function of the neutrino energy is 
shown for the 50(10) 0. 
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Fig. 3 




Figure 3: The survival probability Peri^^e ^e) as function of the neutrino energy is 
shown for the 5*0(10) [0, for sin^(2^eT) = 3.25 x 10~^°. Note that, the corresponding 
value of -E^^, not reported in the figure, results to be much lower than Ea- 
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Fig. 4 




E (MeV) 



Figure 4: The same quantity of Fig. 3 is shown, but with sin^(26'er) = 3.25 x 10 
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Fig. 5 




Figure 5: The energy spectra for are here reported. The sohd hne represents the pure 
initial thermal distribution at the neutrinosphere, whilst the dashed and the dotted lines 
represent the distorted spectrum for 50(10) § with ^\T?{2B^r) = 3.25 x 10"^°, 3.25x 10-^ 
respectively. The dashed-dotted line corresponds to the same quantity for SU{5)susy, 
SO{10)susY, in the DHR scheme @]. 
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Fig. 6 
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Figure 6: The same quantities of Fig. 5 are here shown, but in this case corresponding to 
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